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We introduce the notion of boundary degeneracy of topologically ordered states on a compact 
orientable spatial manifold with boundaries, and emphasize that it provides richer information than 
the bulk degeneracy. Beyond the bulk-edge correspondence, we find the ground state degeneracy 
of fully gapped edge states depends on boundary gapping conditions. We develop a quantitative 
description of different types of boundary gapping conditions by viewing them as different ways 
of non-fractionalized particle condensation on the boundary. Via Chern-Simons theory, this allows 
us to derive the ground state degeneracy formula in terms of boundary gapping conditions, which 
reveals more than the fusion algebra of fractionalized quasiparticles. We apply our results to Toric 
code and Levin- Wen string-net models. By measuring the boundary degeneracy on a cylinder, we 
predict Zk gauge theory and U{l)kXU{l)k non-chiral fractional quantum hall state at even integer k 
can be experimentally distinguished. Our work refines definitions of symmetry protected topological 
order and intrinsic topological order. 

PACS numbers: 



Introduction — Quantum many-body systems exhibit 
surprising new phenomena where topological order and 
the resulting fractionalization are among the central 
themes [1, 2]. One way to define topological order is 
through its ground state degeneracy(GSD) on two spa- 
tial dimensional(2D) higher genus closed Riemann sur- 
face, which is encoded by the fusion rules of fractionalized 
quasiparticles and the genus number [3]. However, on a 
2D compact manifold with boundaries (Fig. 1), there can 
be gapless boundary edge states. For non-chiral topo- 
logical order states, where the numbers of left and right 
moving modes equal, we show there are rules that edge 
states can be fully gapped out. It is the motivation of 
this paper to understand GSD for this type of systems 
where all boundary edge excitations are gapped. In the 
following, we name this degeneracy as boundary GSD, 
to distinguish it from bulk GSD of a gapped phase on 
a closed manifold without boundary. To understand the 
property of boundary GSD is both of theoretical inter- 
ests and of application purpose where lattice model in 
topological quantum computation such as Toric code [4] 
can be put on space with boundaries [5, 6]. 

In this paper, we demonstrate boundary GSD is not 
simply a factorization of the degeneracies of all bound- 
aries. Not only the fusion rules of fractionalized quasipar- 
ticles (anyons) and the manifold topology, but also bound- 
ary gapping conditions are the necessary data. Boundary 
GSD reveals richer information than bulk GSD. More- 
over, gluing edge states of a compact manifold with 
boundaries to form a closed manifold, enables us to 
obtain bulk GSD from boundary GSD. We first intro- 
duce physical notions characterizing this boundary GSD 
and then rigorously derive its general formula by Chern- 
Simons theory [1]. For concrete physical pictures, we re- 
alize specific cases of our result by Z2 Toric code and 
string-net model [7]. Interestingly, by measuring bound- 



ary GSD on a cylinder, our result predicts distinction be- 
tween Zk gauge theory (Z/c Toric code) and U{l)kXU{l)k 
non-chiral fractional quantum hall state at even integer 
/c, despite the two phases cannot be distinguished by bulk 
GSD. By the same idea, we refine definitions of symmet- 
ric protected topological order{SFT) and intrinsic topo- 
logical order. Our prediction can be tested experimen- 
tally. 




FIG. 1: (a) Illustration of fusion rules and total neutral- 
ity, where anyons are transported from one boundary to an- 
other (red arrows), or when they fuse into physical excita- 
tions (blue arrows) (b) A higher genus compact surface with 
boundaries (punctures) . 

Physical Notions — Topological order on a compact 
spatial manifold with boundaries have N branches of 
gapless edge states [1]. Suppose the manifold has to- 
tal T] boundaries, we label each boundary as 9^, with 
1 < a < 77. Let us focus on the case the manifold is home- 
omorphic to a sphere with r] punctures (Fig. 1(a)), we will 
comment on cases with genus or handles(Fig. 1(b)) later. 
If electrons condense on the boundary due to edge states 
scattering, it introduces mass gap to the edge states. A 
set of electrons can condense on the same boundary if 
they do not have relative quantum fluctuation phases 
with each other, thus all condensed electrons are stabi- 
lized in the classical sense. It requires condensed electrons 
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have relative zero Aharonov-Bohm (charge-flux) phase, 
we cah these electrons are null and mutual null[8]. Since 
electrons have discrete elementary charge unit, we label 
them as a dimension- A/" ( dim- A^) lattice Fg, and label con- 
densed electrons as discrete lattice vectors £^°'{g Fg) for 
boundary da- We deflne a complete set of condensed elec- 
trons, labeled as lattice F^°^ , to include all electrons which 
are null and mutual null to each other. Notably there 
are different complete sets of condensed electrons, repre- 
senting different kinds of boundary gapping conditions. 
Assigning a complete set of condensed {non-fractionalized 
bosonic) electrons to a boundary corresponds to assign- 
ing a type of boundary gapping condition. The number 
of types of complete sets equals to the number of types 
of boundary conditions. In principle each boundary can 
assign its own boundary condition independently, this 
assignment is not determined from the bulk information. 
Below we focus on the non-chiral orders, assuming all 
branches of edge states can be fully gapped out. 

Remarkably there exists a set of compatible anyons do 
not produce flux effect to condensed electron charge, so 
their Aharonov-Bohm phases are zero. In other words, 
compatible anyons are mutual null to any elements in the 
complete set of condensed electrons. For boundary 9q,, 
We label compatible anyons as discrete lattice vectors i^^ 
and flnd all such anyons to form a complete set labeled 
as F^^. Note that F^" C F^^. F^« and F^^ have the 
same dimension. If compatible anyons can transport be- 
tween different boundaries of the compact manifold, they 
must follow total neutrality - the net transport of compat- 
ible anyons between boundaries must be balanced by the 
fusion of physical particles in the system(Fig. 1(a)), so 
Sa^gp ^ ^e- Transporting anyons from boundaries to 
boundaries in a fractionalized manner(i.e. not in integral 
electron units), changes topological sectors(i.e. ground 
states) of the system. Given data: Fe,F^°^,F^^, we thus 
derive generally GSD counts the number (Num) of ele- 
ments in a quotient group: 
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GSD = Num[ 



(1) 

Boundary degeneracy of Abelian topological or- 

del — To demonstrate our above notions, let us take 
Abelian topological order as an example, which is be- 
lieved to be fully classifled by K matrix Abelian Chern- 
Simons theory [9]. For a system lives on a 2D compact 
manifold M with ID boundaries dM, edge states of each 
closed boundary (ho meomorphic to S^) are described by 
multiplet chiral bosons[l], with bulk and boundary ac- 
tions: 
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Kjj and Vjj are symmetric integer N x N matrices, 
is the 1-form emergent gauge field's I-th component 
in the multiplet. cos(^^ j ■ ^j) is derived from local 
Hermitian gapping term t/j -\- ip^ = e^^^^^^'^^ + g-*^a,j-^i^ 
where has A'-component with integer coefficients. 

Canonical quantization of K matrix Abelian Chern- 
Simons theory edge states — In order to understand the 
energy spectrum or GSD of the edge theory, we study the 
'quantum' theory, by canonical quantizing the boson field 
^j. Since is the compact phase of a fermion field, its 
bosonization has zero mode ^o/ and winding momentum 
P(f)j, in addition to non-zero modes[10]: 



_i 27r 
^i{x) = (l)oi^KjjP^jj-x- 



(4) 
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The periodic boundary size is L. The conjugate momen- 
tum field of ^i{x) is Ui{x) = j^i^ = i^Kijd:,^j. 
With conjugation relation for zero modes: [^o/7^(/)j] = 
and generalized Kac-Moody algebra for non-zero 
modes: [oLi^n-)^j,m\ = ^Kjjdn^-ra- We thus have canon- 
ical quantized fields: [<l>/(xi), I\.j{x2)] = \i8ij8{xi — X2). 

Gapping Rules, Hamiltonian, and Hilbert Space [11] — 
Let us determine the properties of ^a as condensed elec- 
trons and the set of allowed boundary gapping lattice 
= {L}, ci labels the a-th vector in F^. With- 
out any symmetry constraint, then any gapping term 
is allowed if it is [11]: (1) Null and mutual null [8]: 
V£^,£^gF^, ^6,j = 0' This implies seff statis- 

tics and mutual statistics are bosonic, and the excitation 
is local. Localized fields are not eliminated by self or 
mutual quantum fiuctuations, so the condensation sur- 
vives in the classical sense. (2) 'Physical' excitation: 
£^ G Fg = {^j cjKij I cj e Z}, is an excitation of 
electron degree of freedom since it lives on the 'physical' 
boundary. (3) Completeness of F^, defined by: \/£^ G Fg, 
if i^K-^i^ = and ic^-^L = with G F^, then 
£^ G F^ must be true. It turns out rules (1)(2)(3) are not 
guaranteed to fully gap out edge states, we add an extra 
rule (4): The system is non-chiral. The signature of K 
(= Num of positive eigenvalues — Num of negative eigen- 
values) must be zero, i.e. A^ is even. Additionally, the 
dimension of F^ must be N/2 and det K\ G N. Phys- 
ically, rule (4) excludes violating examples such as odd 
rank(=rk) K matrix with central charge c — Cy^ 0, which 
universally has gapless chiral states. For instance, the 
dim-1 boundary gapping lattice: {n{A, B,C) | n G Z} of 
^3x3 = diag{l, 1, -1}, with + 5^ - = 0, satisfies 
(1)(2)(3), but cannot fully gap out edge states. 
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Since any linear combinations of G Fg still sat- 
isfy (1)(2)(3)(4), we can regard F^ as an infinite dis- 
crete lattice group. As we solve exactly the number of 
types of boundary gapping conditions A/"^, at rk(i^) = 
2, AT^ = 2. However, for rk(K) > 4, = oo, 

which is rather surprising. For example, as i^4x4 = 
diag{l, 1,-1,-1}, one can find dim-2 boundary gapping 
lattices, 5, C, 0), m(0, C, 5, A) | n, m G Z}, where 

different sets of + 5^ — = give different lattices. 
There are infinite solutions of them[ll]. 

To determine the mass gap of boundary states and 
show the gap is finite in the large system size limit 
1/ ^ oo, we will take large g coupling limit of Hamilto- 
nian: -ga dx cos{£^ j-^i) ^gaii^j - ^i)'^L. Solve 
the quadratic Hamiltonian exactly, we find the spectrum 
E{n) = + + (^)(52, generically indicating 

edge states has finite gap A independent of boundary size 
L. To count boundary GSD of this gapped system, we 
need to include the full cosine term for the lowest energy 
states - zero and winding modes: 

cos(^^,j • ^i) ^ cos(^^^^ • (^0/ + K^JPc^j -j-x)) (5) 



The kinetic term Hkin = ■^^^/J^i^{^7z2^0zi^0^2 has 
0{1/L) so can be neglected in L ^ oo, the cosine po- 
tential Eq. (5) dominates, use [^ai^oij i'-^r^j-^(/>j] ~ 
i^T'^j^a I'^iJ = by rule (1), we can safely expand 
the cosine term, then integrate L\ 



over all r] boundaries must be on the Fg lattice. 



ga / dx Eq.(5) gaL cos{£^ j 
Jo 



(6) 

^0 is periodic, so forms a discrete lattice. By rule (2), 
so cos(£^ / • 00/) are hopping terms along condensed elec- 
tron vector j in sublattice of F^ in lattice. P^^ rep- 
resents compatible anyons £qp which is mutual null to con- 
densed electrons £ by £K-^Pf = £K-^£qp = 0. By rule 
(1), £qp parallels along some £ vector. However, £qp lives 
on quasiparticle lattice, i.e. unit integer lattice of P^ lat- 
tice. So £qp is parametrized by |gcd(£ )\ ^a,j^ ^^^^ greatest 
common divisor | gcd(4)| = gcd(|4,i|, |4,2|, • • • , |4,Ar|). 

For boundary da, a complete set of condensed elec- 
trons forms boundary gapping lattice F^°^ = {J^a ^a'^^a'^i I 
/f" G Z}. A complete set of compatible anyons forms the 
Hilbert space of winding modes on P^ lattice: 



£' 



a, I 



?cd(£^)| 



fa" e Z} (7) 



or simply anyon hopping lattice. Anyon fusion rules and 
total neutrality essentially imply physical charge excita- 
tions can fuse into multiple anyon charges. The rules 
constraint the direct sum[12] of anyon hopping lattice F^^ 



?cd(£^)| 



Vjf" G Z, 3 cj G 



^ |gcd(4-)| Y 

GSD counts topological sectors distinct by fractionalized 
anyons transport between boundaries. With 0^ F^°^ is 
a normal subgroup of PgpQe^ and given input data: K 
and F^" (which determines F^^), we derive GSD is the 
number of elements in a quotient finite Abelian group: 



GSD = Num [ 



(9) 



analog to Eq. (1). Interestingly Eq. (9) works for 
both closed manifolds or compact manifolds with 
boundaries[13]. By gluing boundaries of compact man- 
ifold, enlarging K matrix of glued edges states to 
^2 AT X 2 AT, creating N scattering channels to fully gap 
out edge states. For a genus g Riemann surface with 
T]^ punctures(Fig. 1(b)), we start with a number of g 
cylinders drilled with extra punctures [11], use Eq. (9), 
remove glued boundaries part which contributes a factor 
< I deti^l^, and redefine particle hopping lattices Pgpf|e 
and 0^/ F^«' only for unglued boundaries {I < < 77'), 



GSD< |deti^|^-Num[ 



^qpf\e 



(10) 



For genus g Riemann surface (77' = 0), Eq. (10) becomes 
GSD < |deti^|^. The inequalities are due to different 
choices of gapping conditions for glued boundaries [13]. 

Apply our algorithm to generic ^^2x2 matrix case [11], 
to fully gap out edge states require detK = —k'^ with 
integer k. Take a cylinder (a sphere with 2 punctures) as 
an example, Eq. (9) shows GSD = det K\ = k when 
boundary conditions on two edges are the same F^^ = 
However, GSD < k when F^^ ^ Specifically, 
take = ( ^ ) ^diag,/c = ( -/c ) ' we obtain Kz^ 
has GSD = 1, while i^diag,/c has GSD = 1 for odd k but 
GSD = 2 for even k. See Table F 







-^diag,fc 


Boundary 




1 


1 (/c G odd) or 2 (/c G even) 


GSD 




k 


k 


Bulk GSD 







TABLE L Boundary GSD and bulk GSD for Kz^ and Kaie.g,k 

This shows a new surprise, we predict a distinction 
between two types of order Kz^i^k gauge theory) and 
Kdiag,k{U (l)/c X /7 {l)k non-chiral fractional quantum hall 
state) at even integer k by simply measuring their bound- 
ary GSD on a cylinder. 
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Mutual Chern- Simons, gauge theory , Toric 
code and String-net model — We now take Kzf, ex- 
ample to demonstrate our understanding of two types of 
GSD on a cylinder in physical pictures. 



Chern-Simons Toric Code String-net 




FIG. 2: (a) Same boundary conditions on two ends of a cylin- 
der allow a pair of cycles [c^], [cz] of a qubit, thus GSD = 
2. Different boundary conditions do not, thus GSD = 1. 
(b)Same boundary conditions allow z- or x-strings connect 
two boundaries. Different boundary conditions do not. 

When /c = 2, it realizes Z2 Toric code [14] with 
Hamiltonian[4] Hq = — Ay — Bp. There are two 
types of boundaries [5] on a cylinder (Fig. 2(a)): x (R: 
rough) boundary where z-string charge e condenses and 
z (S: smooth) boundary where x- string charge m con- 
denses. We can determine GSD by counting the de- 
gree of freedom of the code subspace: Num of qubits 
— Num of independent stabilizers. For F^^ = F^^^ 
have the same number of qubits and stabilizers, with 
one extra constraint Yl^n Bp = 1 for two x-boundaries 
(ria^^ = 1 for two 2:-boundaries) . This leaves 1 free 
qubit, thus GSD = 2. For F^^ 7^ F^^ stih the same num- 
ber of qubits and stabilizers, but has no extra constraint. 
This leaves no free qubits, thus GSD = 1. We can also 
count independent logical operators(Fig. 2(a)) in homol- 
ogy class, with string-net (Fig. 2(b)) picture in mind - 
there are two cycles [c^^J, [c^J winding around the com- 
pact direction of a cylinder. If both are x-boundaries, we 
only have z-string connects two edges: cycle [cz^]- If both 
are z-boundaries, we only have x-string (dual string) con- 
nects two edges: cycle [0x2]- Cycles of either case define 
the algebra cF^^cFy of a qubit, so GSD = 2. For both 
x-boundaries (z-boundaries), one ground state has even 
number of strings (dual strings), the other ground state 
has odd number of strings (dual strings), connecting two 
edges. If boundaries are different, no cycle is allowed in 
the non-compact direction, no string and no dual string 
can connect two edges, so GSD = 1. Generally, for a 
level k doubled model, GSD = k'^^M^iiM,'^^)] ^ k^iiM) ^ 
k to the power of the 1st Betti number [13, 15]. 
Definitions of Topological Order — Now let us ask a 
fundamental question: what is topological order? We re- 
alize the original definition of degenerated ground states 



on a higher genus Riemann surface can be transplanted to 
degenerated ground states on a cylinder with two bound- 
aries. For a non-chiral fully-gapped-boundary state, we 
define: the state is intrinsic topological order, if it has 
degenerated ground states ( at least for certain boundary 
gapping conditions) on a cylinder. The at /ea^t statement 
is due to GSD < det K|, only the same boundary con- 
ditions on two sides gives GSD = ^\deiK\. Similarly, 
the state is trivial order or SPT, if it has a unique ground 
state on a cylinder for any boundary gapping condition. 
Conclusion — In summary, by K^xN matrix Abelian 
Chern- Simon theory, we derive boundary fully gapping 
conditions for Abelian topological order states, its low 
energy Hamiltonian and Hilbert space, a general GSD 
formula, and the number of types of boundary gap- 
ping conditions JVg . The fully gapped boundaries re- 
quires N = even and non-chiial^doubled) Chern-Simons 
theory, reflecting gapped boundary Quantum Doubled 
model of Toric code and string- net [6, 16]. Our JVg for- 
mula proves the two boundary types conjecture of Toric 
code [5, 15, 16], Zk gauge theory, and a more general ^2x2 
Chern- Simon theory. We show gapping edge states to 
count boundary GSD is related to bulk GSD. However, 
we find there are more types of boundary GSD (than a 
unique bulk GSD) depending on types of boundary gap- 
ping conditions. A remarkable example is rank- 2 Kzj, 
and i^diag,/c with k = even[ll], where we predict exper- 
imentally GSD = A: or 1 for Kz,^ and GSD = A: or 2 for 
^diag,/c on a cylinder, though GSD of both phases on a 
closed genus g surface are indistinguishable (= k'^^). This 
example is especially surprising because both phases Kzj, 
and i^diag,/c have the same (Z/^)^ fusion algebra. This 
means fusion algebra alone cannot determine boundary 
GSD [11]. In the category language [6], the model of uni- 
tary fusion category C shows that (in Table II) there can 
be many different C, V types realizing the same monoidal 
center Z{C) = Z(V). 



Physics 


Category 


Bulk excitation 


objects in unitary modular category 


(anyons) 


Z{C) (monoidal center of C) 


Boundary type 


the set of equivalent classes {C, P, . . . } 
of unitary fusion category 



TABLE II: Dictionary between physics and category 



Finally, our definitions of topological orders not only 
deepen understanding of topological GSD, but also ease 
the experimental platform with only cylinder topology 
instead of higher genus surfaces. To open future re- 
search avenues, it will be interesting to realize boundary 
GSD and Mg without using K matrix, which is restricted 
only to Abelian topological order, also to study them on 
higher dimensional space. Physical notions which we in- 
troduced should still hold universally. Other than fusion 
rule and total neutrality^ whether braiding rule can enter 
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into GSD formula[17]? Additionally, engineer boundaries 
to line-like defects may synthesize projective non-Abelian 
statistics [18]. All these shall inspire generalizing Eq. (9) 
to boundary GSD of non-Abelian topological order. 

It will be illuminating to have more predictions based 
on our theory, as well as experimental realizations of 
boundary types. One approach is flux insertion through 
a cylinder, where the adiabatic flux change induces 
anyon transport from one boundary to the other by 
A^B/{h/e) = AP0[11]. It win be interesting to see 
how the same type of boundary conditions allow this ef- 
fect (unit flux changes topological sectors, with total sec- 
tors as GSD = det K\)^ while different type of bound- 
ary conditions restrain this effect 'dynamically' (GSD < 
y^l det K\). To detect this dynamical effect can guide 
experiments to distinguish boundary types. 
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APPENDIX 

In Appendix, we explain detailed physical meanings of 
boundary gapping rules in Sec. A, and demonstrate our 
algorithm and GSD formula Eq. (9) for a generic rank-2 
K matrix in Sec.B. We also give an example why fusion 
algebra alone does not provide enough information to 
determine boundary GSD from the bulk-edge correspon- 
dence viewpoint. In Sec.C, we comment more about the 
number of boundary types JVg . In Sec.D, we use gluing 
techniques to derive Eq. (10) for a compact manifold with 
genus. Lastly, we explain how flux insertion experiment 
may help to test boundary types in Sec.E. 



A. BOUNDARY GAPPING RULES 

We argue boundary gapping rules are: 

(i) Bosonic: V^a ^ L^, jKjji^ j G even integer^ 
which means mutual statistics has multiple 27r phase. 

(ii) Local: V^^,^^ G F^, ia,i^7j^b,j ^ integer. Winding 
one excitation around another yields multiple 27r phase. 

(iii) Localizing field at classical value without being 



eliminated by self or mutual quantum fiuct nation: 
V^^,^5 G F^, £^ jKjj£^ J — ^^ so condensation 
survives in the classical sense. 

(iv) must be excitations of electron degree of freedom 
since it lives on the 'physical' boundary, so G Fg elec- 
tron lattice, where Fg = {Xlj cjKjj \ cj G Z} . This rule 
imposes integer charge qjKjjl^ j in the bulk, and in- 
teger charge Qi = ^d^^^idx = KjjP^j = K'ji^ j 
for each branch on the boundary. Here qj is the charge 
vector coupling to an external field of gauge or global 
symmetry, by adding A^qjJj to the Shuik-, which causes 
(lI^^l^x^I in the Sq- 

(v) Completeness: we define F^ is a complete set by: 
V4 G Fe, if ic^-^^c = and ^i^"^^ = with 
Wi^ G F^, then £^ G F^ must be true, otherwise F^ is not 
complete. The completeness rule is due to no symmetry 
protection forbids any possible local gapping terms. 

We can summarize all above rules of F^ as (1) Null and 
mutual nulh Vr,^^GF^, ^a,i^7jh,J = (2) 'Phys- 
ical' excitation: £a ^ Fg, (3) Completeness of F^. 
(1)(2)(3) are not guaranteed to fully gap out edge states, 
we add an extra rule: (4) N G even. Dimension of F^ 
must be which requires ^^^det^^ G N. 



B. EXACT RESULTS OF ^2x2 CHERN- SIMONS 
THEORY 

Here we work through a rank-2 K matrix Chern- 
Simons theory example, to demonstrate our generic algo- 
rithm in the main text. We derive its low energy Hamil- 
tonian, Hilbert space, boundary GSD formula, and the 
number of types of boundary gapping conditions A/"^. For 

K2.2 = itl ll)^ ill ). in order to fully gap 

out edge states, we find that deti^ = — p^, (a) the edge 
states need to be non-chiral (i^2x2 has equal number of 
positive and negative eigenvalues, so deti^ < and (b) 
|deti^^| needs to be an integer p square. Two indepen- 
dent sets of allowed gapping lattices F^ = {n£^ j\n G Z}, 
F^' = {n'iayW ^ ^} satisfy gapping rules (1)(2)(3)(4), 
so JV^ = 2 at Tk{K) 2, with 



n p 



\gcd{ki,k3 
p 

|gcd(A:3 +p, k2)\ 



{ki.ks^p) (11) 
{ks^p,k2) (12) 



|gcd(A:,/)| stands for the greatest common divisor in 
its absolute value. If k (or /) is zero, we define 
|gcd(/c,/)| is the other value |/| (or |A:|). n^n^ e Z are 
allowed if no other symmetry constrains its values. 

Now take specific topology, a disk(a sphere with 1 
puncture) and a cylinder(a sphere with 2 punctures) as 
examples of manifolds with boundary, for i^2x2 their 
Hilbert space of edge states are[l]: 

nd^sk = n^K^M ^ ^p,, ^ ^p,, (13) 

(if i^^2x2 is diagonal, then ^km ^^km)- 

a^sk^'^Tr®'H,i (15) 
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1-Lkm stands for Hilbert space of nonzero Fourier modes 
part with Kac-Moody algebra. We label zero and wind- 
ing modes' Hilbert space by winding mode P^, which 
can be regarded as a discrete lattice because of ^/(x) 
periodicity. Because the bulk cylinder provides channels 
connecting edges states of two boundaries, so fractional- 
ized quasiparticles (here Abelian anyons) can be trans- 
ported from one edge to the other. 1-Lgi contains frac- 
tional sectors Ija, Jb), the 1st branch ja runs between 
the top(P^7) and the bottom(P^f ^^^), the 2nd branch 
jB runs between the top(P^7) and the bottom(P^f ^^^). 

Let us explicitly show that edge states of any bound- 
ary has a finite gap from ground states at large system 
size L and large coupling g. Without losing general- 
ity, take V = {v\vl) and a gapping term {^^ai^^a2) — 



|gcd(fef,fc3+p)| (^i^ ^3+p) e r^, and diagonalize the Hamil- 
tonian, 

H^ij^ dxVud:,^id,^j) + ]^g{l,j-^ifL (16) 
we find eigenvalues: 



^i,2(n) = ^A2 + (^)2^,±(^)5, (17) 

where the finite mass gap is independent of L: 

^2Tig{ki(ki - k2)vi + 2(/c3 ^ p)(k2,vi - kiV2)) 
|gcd(A;i,/c3 +p)| 

^1 = (^1 - k2Yvl^A{k^vi - kiV2){ksVi - /c2^2) and 62 = 
vi{ki + k2) - ^2(2/^3). 

To count boundary GSD, for generic i^2x2 Abelian 
topological order on a disk, take £^ G F^ in Eq. (11) with- 
out losing generality (same argument for F^ )[2], we have 



^01 =h^aA^jA\ 



- |gcd(^^)| ' — A^a,2+JA |gcd(£^)| 

total anyon charge for each branch needs to conserve, but 
a single boundary of a disk has no other boundaries to 
locate transported anyons. This implies: Ja = 0, there 
is no different topological sector induced by transporting 
anyons, thus GSD = 1. 

On the other hand, if the topology is replaced to a 
cylinder with the top di and the bottom 82 boundaries 
in Fig. 3, when the gapping terms from boundary gapping 
lattice F^ are chosen to be the same, the Hilbert space 
on P^ lattice is: 



. The 



Pet's = hini+JB-, 



?cd(4)r 
?cd(4)| 



, P^^ = /2^a,2 + JB I 



?cd(4)r 

^a,2 



?cd(4)| 



(19) 



Anyon fusion rule and charge conservation for each 
branch constrains (P^^ + P^p^^P^j)^ + P^^) G Fg electron 
lattice. With | gcd(£a)| = it implies Ja = —jsi'modp). 
^ jAi'^od p) < p has p different topological sectors 
induced by different Ja- When AJa/p G Z, it transports 
electrons, so it brings back to the same sector. Count 
the number of distinct sectors, i.e. ground states, we 
find GSD = p. 

If gapping terms on two boundaries of a cylinder are 
chosen to be different: F^ for 9i, F^ for 82, we revise 
the second line of Eq. (19) to 



^2^,1 



-Jb 



a, 2 



gcd(C)|'^*^ '2'°-2+^^|gcd(C)| 

(20) 

rules and conservation implies: 

^|gcd(£^)K^^ |gcd(£;)KB' |gcd(£^)K^^ IgcdC^'^)^^^ t J-e- 

This constraint gives a surprise. For exam- 
pie, K. - ('^ 



Anyon fusion 

( ' 



(^g), GSD = 1. However, when 
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Kdmg,p = (o -p)^ GSD = 1 for p G odd, but GSD = 2 
for p e even. This provides a new approach that one can 
distinguish two types of orders and i^diag when p is 
even by measuring their boundary GSD. 

We ihustrate this result in an intuitive way in Fig. 3. 
When boundary types are the same on two sides of the 
cyhnder, Fig. 3(a) is enough to explain GSD = p , where 
fractionalized anyons transport from the bottom to the 
top. For case, say F^^ = F^^ ^ = {n(j9,0)}, qpi 
with ia = 0) for < Ja < P - I- For i^diag,p case, 
say F^^ = F^2 = F^ = {n{p,p)}, qpi with 4 = 1) 
for < ^ P — 1- This accounts all p sectors. 

When boundary types are different. Fig. 3(a) is not al- 
lowed for fractionalized anyons transport. Fig. 3(b) is 
crucial to account the second ground state of i^diag,p at 
p G even. Let us take F^^ = F^ = {n{p^p)} and F^^ = 
F^ = {n{p^ —p)}^ where e^ represents £ = (p, 0), while 
qpi with 4 = (p/2,p/2) and qp2 with £a = (p/2, -p/2) 
at p G even are allowed fractionalized anyons (with in- 
teger unit of anyon charge). This process switches the 
ground state to a different sector, so GSD = 2. However, 
fractionalized anyons transport in Fig. 3(b) is not allowed 
for with different boundary types on two sides of the 
cylinder, which results in GSD = 1. 




FIG. 3: (a)Anyon(qpi) is transported from the bottom to 
the top of the cylinder. (b)Physical excitations e^ split into a 
pair of anyons(qpi to the bottom, qp2 to the top). Note that 
(b) is crucial to explain GSD = 2 for i^diag,p at p G even with 
different boundary types on two sides of the cylinder. 

The result remarks that only fusion algebra(both have 
doubled fusion algebra (^p)^) is not sufficient enough 
to determine boundary GSD. Let us show this by 
and i^diag,2 examples from the bulk-edge correspondence 
viewpoint. Both orders have the doubled fusion al- 
gebra, (^2)^, with bulk excitation such as two anyon 
types 4 = (0,1), 4 = (1,0). For Kz^, two types 
of boundary conditions F^ = {n2(0, 1) | n G Z} and 
F^ = {n2(l,0) I n G Z} correspond to two of 4,4, thus 
this is one to one correspondence. However, for i^diag,2, 
two boundary conditions F^ = {n2(l, 1) | n G Z} and 
F^ = {n2(l,-l) I n G Z} correspond to one bulk ex- 
citation i = {l{mod2), l{mod2)) = {l{mod2), -l{mod2)) 
in fusion algebra. It is two (on the boundary) to one (in 
the bulk) correspondence. Hence i^diag,2 shows an ex- 
ample where anyon fusion algebra from the bulk cannot 
distinguish different boundary types. 



C. NUMBER OF TYPES OF BOUNDARY 
GAPPING CONDITIONS 

For the number of types of boundary gapping condi- 
tions J\f^, at Tk{K) = 2, we showed = 2, which proves 
the conjecture [3] nicely, rk(i^) > 4, Mg is infinite^ is sur- 
prising. We may wonder whether extra rules are possible 
to restrict J\fg . Let us assume |deti^| = 1, the canon- 
ical form of this unimodular indefinite symmetric inte- 
gral KnxN matrix exists [4]. For the odd matrix(where 
quadratic form has some odd integer coefficient, so with 
fermionic statistics), the canonical form is composed by 
N/2 blocks of ( -^1 ) along the diagonal blocks of KnxN- 
For the even matrix (where quadratic form has only even 
integer coefficient, so with only bosonic statistics), the 
canonical form is composed by blocks of ( 1 J ) a set 
of all positive(or negative) coefficients £^8 lattices, Kes^ 
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(21) 



along the diagonal blocks of KnxN- A positive defi- 
nite Kes with eight chiral bosons cannot be gapped out. 
Thus, in order to have non-chiral states, the even matrix 
canonical form must be composed by N/2 blocks of (10) 
Now it is useful to rethink about the number of boundary 
types J^g in this canonical form when | det K\ = 1. We 
had claimed under rule (1)(2)(3)(4) when rk(i^) > 4, 
^fg =00. If we add an extra gapping rule (5): Each 
gapping term is formed by only two branches of chiral- 
antichiral edge states scattering. There are no three or 
more multi-channel scatterings. With (5), the bound- 
ary gapping lattices F^ are formed by a set of £ terms 
of two branches of edge states with opposite chirality. 
This modifies our result to J\f^ = {N/2)\2^/^. Its {N/2)\ 
factor is due to each chiral state can find another state 
from antichiral ones to form a pair. The relative sign of 
their bosonized phase is ±, gives an extra 2^/^ factor. 
We comment that (5) seems to be artificial without deep 
reasoning, and the formula only works for |deti^| = 1 
with the canonical form. However, one can postulate 
some extra rule like (5) can restrict Afg . 

D. SURGERY TO GLUE CYLINDERS TO FORM 
A GENUS g RIEMANN SURFACE WITH 
PUNCTURES 

Here we show how to glue the boundaries of punctured 
cylinders to form a genus g Riemann surface with punc- 
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tures, and determine its GSD of the topological order. 

For a genus g Riemann surface with rj^ punc- 
tures(Fig. 1(b) and Fig. 4(c)), we start from Fig. 4(a), a 
number of g cylinders drilled with total puncture num- 
ber T] = T]' -\- 2g -\- 2{g — 1), where 2g count two punc- 
tures on top and bottom for each cylinder(hi^T and h^^^, 

1 < ^ < ^) , drill an extra puncture on both the 1st 
(hi,L) and the last ^th cylinder (h^_i^i^), and drill two 
extra punctures (hj_i^i^ and ^j^l) for j-th cylinder for 

2 < j < g — 1. There are thus 2{g — 1) extra punc- 
tures. Glue the boundaries of hj^L and hj^ji together 
for 1 < j < ^ — 1, and glue h^^^ and h^^^ together for 
1 < ^ < ^, results in Fig. 4(c). Use Eq. (9), remove glued 
boundaries part(l < a < 2g-\-2{g — l)) which contributes 



a factor < \detK\^^ and redefine particle hopping lat- 
tices Lqp pi e and 0^, r^«' only for unglued boundaries 
(1 < a' < 7^'), we derive 

GSD < I det K\^ • Num [ ] (22) 

For a genus g Riemann surface (r^^ = 0), this becomes 
GSD < I det K\9^ where rk(i^) = AT for a closed manifold 
case is relaxed to any natural number N, which works for 
both odd and even number of branches. The inequalities 
here are due to different choices of gapping conditions for 
glued boundaries. 



hl,T ^ h2,T , , .^ hg-l,T hg,T 




FIG. 4: Glue punctured cylinders to form a genus g Riemann surface with r]' punctures. Start from (a), firstly identify left and 
right > arrows of each square to form number g of punctured cylinders. Then glue h^,!, and hj^R (red dotted circles) together 
for 1 < j < ^ — 1, and glue hi^r and hi,s (blue arrows) together for 1 < i < ^, which yields (b), equivalently as a genus g 
Riemann surface (c). The extra r/' punctures are indicated here as a shaded blue puncture in the left most handle. 



E. FLUX INSERTION ARGUMENT AND 
EXPERIMENTAL TEST ON BOUNDARY TYPES 

Consider an artificial-designed or external gauge field 
(such as electromagnetic field) coupled to topologically 
ordered states by a charge vector qj. An adiabatic fiux 
insertion A^b inside the cylinder induces E^^ causing 
a perpendicular current Jy fiows to the boundary, from 
the bulk term = -qi^Kjj ^e^^^PdyAp, so qiA^B = 
-qi J dt jE-dt= -^Kijh J Jy^jdtdx = -^Kij^Qj. 
Qj is the charge condensed on the edge of the cylin- 
der. On the other hand, the edge state dynamics affects 
winding modes by Q/ = J Jq jdx = — J ^dx^idx = 
-eKjjP^^j, so 

qiA^uK-) = AP^j (23) 



For the same types of boundaries: F^^ = F^^^ there are 
^y\detK\ sectors by AP^(mod vT^et^), the ^J\dGiK\ 
units of fiux bring the state back to the original sector. 
For F^i ^ r^2^ we had shown GSD < ^\d^iK\ (such 
as GSD = 1). This motivates an interesting question 
if one inserts fiux into the cylinder, what dynamical ef- 
fect, which repulses anyons transporting from one edge 
to the other, will be detected. The detection of this effect 
may guide experiments to distinguish different types of 
boundary gapping conditions. 



[1] X. G. Wen, Phys. Rev. B 41, 12838 (1990). 

[2] Notice if we take |n| > 2 or \n'\ > 2 larger-size hopping 
term for niaj G ,'n^^a,i ^ F^ , the lattice can break 
into more sublattices of different sectors: the calculation 
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on a disk is altered to GSD = \n\. However, generically [3] S. B. Bravyi, A. Y. Kitaev, arXiv: quant-ph/9811052. 
no symmetry forbids n = = 1 least-size hopping terms, [4] littp://matlioverflow. net /questions/97448 
which dominant cosine potential can minimize the energy. 



